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1.  Introduction 

Consider  a  square  with  area  N  in  the  plane  that  contains  N  grid  points  arranged 
in  a  regularly  spaced  y/N  x  y/N  array  and  N  random  points  located  independently  and 
randomly  according  to  the  uniform  distribution  on  the  square.  For  example,  see  Figure  1. 


Figure  1:  A  square  with  area  N  containing  a  regularly  spaced  y/N  x  y/N  grid  (denoted 
with  black  dots)  and  N  uniformly  distributed  random  points  (denoted  with  x’s). 

For  any  particular  set  of  N  random  points  P,  let  L(P)  denote  the  minimum  length 
such  that  there  exists  a  perfect  matching  of  the  (random)  points  in  P  to  the  grid  points  in 
the  square  for  which  the  distance  between  every  pair  of  matched  points  is  at  most  L(P). 
In  other  words,  L(P)  is  the  minimum  over  ail  perfect  matchings  of  the  maximum  distance 
between  any  pair  of  matched  points,  or  more  simply  the  minimal  matching  length  for 
P.  As  an  example,  Figure  2  illustrates  two  matchings,  one  that  achieves  the  minimax 
matching  length  and  one  that  does  not. 


Figure  2:  Two  possible  matchings  for  a  set  of  random  points,  one  which  achieves  the 
minimax  matching  length  ( b )  and  one  which  does  not  (a). 


Algorithmically  speaking,  the  problem  of  computing  the  minimax  matching  length  for 
any  set  P  is  relatively  straightforward,  and  is  not  of  concern  in  this  paper.  Nor  are  we 
interested  in  the  worst  case  value  of  L(P)  over  all  P,  which  is  trivially  0(\/N).  Rather, 
we  are  interested  in  the  expected  value  of  L(P)  for  random  P. 

Initially,  one  might  hope  that  the  expected  value  of  L(P)  is  a  constant  independent  of 
N  since  every  random  point  is  within  distance  y/2/2  of  a  grid  point.  With  high  probability, 
however,  the  same  is  not  true  for  every  grid  point.  In  particular,  it  is  not  difficult  to  show 
that  with  probability  exceeding  1  —  jf  there  is  a  circular  region  with  ©(log#)  area  in  the 
square  that  does  not  contain  any  random  points  at  all,  and  hence  L(P)  >  fi(v/I°g  N)  with 
high  probability.  For  example,  see  Figure  3. 


Circle  of  uu 
f (log  N)  not 
containing  nny 
random  point* 


Figure  3:  Collection  of  N  random  points  P  for  which  there  is  a  circle  of  area  ©(log  N) 
not  containing  any  random  points.  The  grid  point  (not  shown  here)  nearest  to  the  center 
of  this  circle  is  fl(\/log N)  away  from  every  random  point  and  hence,  L(P)  >  ft(y/\ogN) 
for  this  example. 


Although  the  minimal  grid  matching  problem  (that  of  determining  the  expected  value 
of  L(P)  for  random  P  ),  has  not  to  our  knowledge  been  directly  considered  before,  Leighton 
and  Leiserson  [13]  and  Ajtai,  Koml6s  and  Tusn&dy  [1]  considered  similar  problems,  and 
developed  probabilistic  divide-and-conquer  techniques  that  show  L(P )  <  O(logiV)  with 
high  probability.  The  resulting  G(y/Tog77)  gap  in  the  bounds  for  the  expected  value  of 
L(P)  remained  for  some  time  until  Shor  [21,  22]  proved  that  L(P)  >  n(log5^4  N)  with  very 
high  probability *  (i.e.,  with  probability  exceeding  1  —  1  /Na  where  a  =  n(v/log  N)).  Shor's 
result  was  substantially  more  difficult  than  the  rather  elementary  fi(\/log  N)  bound,  and 
suggested  that  the  minimax  grid  matching  problem  possessed  a  deeper  structure  than  one 
might  initially  realize. 


i 


In  this  paper,  we  complete  the  asymptotic  analysis  of  the  minimax  grid  matching 
problem  by  proving  that  L(P)  <  0(loga^4N)  with  very  high  probability.  Hence  L(P)  = 
0(log3/4  N)  with  very  high  probability,  and  the  expected  value  of  L{P)  is  0(logs^4  JV). 
Aside  from  being  an  interesting  combinatorial  problem,  the  improved  upper  bound  for 

*  Henceforth,  we  will  reserve  the  phrase  *with  high  probability*  to  mean  ‘with  probability  exceeding 
1  -  1  /JVU  for  any  constant  a  >  1,*  and  the  phrase  “with  very  high  probability*  to  mean  “with  probability 
exceeding  1  —  \/Na  where  a  =  n(>/log  N). 
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typical  values  of  L(P)  has  important  implications  for  the  average  case  behavior  of  a  wide 
variety  of  algorithms.  For  example,  as  a  direct  consequence  of  this  result,  we  now  know  that 
the  1-dimensional  Best  Fit  bin  packing  algorithm  wastes  0(\/77logs^4  N)  space  with  very 
high  probability  when  packing  N  items  with  sizes  determined  by  identical  and  independent 
distributions  that  are  symmetric  about  1/2  in  [0,1].  Similar  consequences  hold  for  the 
Karp-Luby-Marchetti-Spaccamela  2-dimensional  bin  packing  algorithm  and  the  Coffman- 
Leighton  Best  Fit  Aligned  dynamic  allocation  algorithm.  Our  result  also  has  implications 
for  the  maximum  up-right  matching  problem,  a  long  open  discrepancy  problem,  wafer-scale 
integration  of  systolic  arrays,  and  potentially  for  testing  pseudorandom  number  generators. 

Establishing  the  connection  between  the  minimax  grid  matching  problem  and  the  wide 
variety  of  seemingly  unrelated  problems  just  mentioned  has,  of  course,  required  a  great  deal 
of  work  by  many  researchers.  Hence,  it  is  not  possible  to  provide  detailed  explanations  for 
all  of  the  connections  in  this  paper.  Instead,  we  include  only  the  relevant  definitions  and 
statements  of  results  for  each  connection,  referring  the  interested  reader  to  other  sources 
for  details. 

The  remainder  of  the  paper  is  divided  into  six  sections.  In  Section  2,  we  describe  some 
related  combinatorial  problems.  In  particular,  we  discuss  the  maximum  up-right  matching 
problem  and  several  2-dimensional  discrepancy  problems.  In  Section  3,  we  discuss  the 
application  of  the  minimax  grid  matching  problem  to  bin  packing,  dynamic  allocation, 
wafer-scale  integration  and  testing  pseudorandom  number  generators.  The  proof  of  the 
0(log3/4  N)  upper  bound  for  minimax  grid  matching  is  contained  in  Section  4.  We  conclude 
with  some  remarks,  acknowledgements  and  references  in  Sections  5-7. 

2.  Related  Combinatorial  Problems 

In  this  section,  we  discuss  several  combinatorial  problems  that  are  closely  related  to  the 
minimax  grid  matching  problem.  These  problems  are  interesting  in  their  own  right,  but  are 
included  here  mostly  for  their  role  in  the  applications  described  in  Section  3  and  the  proof 
contained  in  Section  4.  We  commence  with  the  maximum  up-right  matching  problem  in 
Section  2.1  and  finish  with  a  discussion  of  two-dimensional  discrepancy  problems  in  Section 
2.2. 

2.1.  The  Maximum  Up-Right  Matching  Problem 

Consider  a  square  with  area  N  in  the  plane  that  contains  N  random  pluses  and  N 
random  minuses  located  independently  and  randomly  according  to  the  uniform  distribution 
on  the  square.  Given  a  set  of  N  pluses  and  N  minuses  P*,  an  up-right  matching  on  P± 
is  a  one-to-one  matching  of  pluses  to  minuses  such  that  every  plus  is  either  unmatched  or 
is  matched  to  a  single  minus  that  lies  above  and  to  the  right  of  the  plus.  A  maximum  up¬ 
right  matching  for  P±  is  an  up-right  matching  which  minimizes  the  number  of  unmatched 
points.  For  example,  Figure  4  illustrates  two  up-right  matchings.  One  is  maximum,  the 
other  is  not. 

Let  U(P±)  denote  the  number  of  unmatched  pluses  in  a  maximum  up-right  matching 
for  P±.  For  example,  U(P±)  =  1  for  the  collection  of  pluses  and  minuses  shown  in  Figure 
4.  As  with  minimax  grid  matching,  the  algorithmic  problems  of  determining  U(P±)  and 
finding  a  maximum  up-right  matching  for  a  particular  collection  of  pluses  and  minuses  are 


Figure  4:  Two  up-right  matchings.  Only  the  matching  in  (b)  is  maximum. 


elementary  and  not  of  interest  here.  Rather,  we  are  interested  in  determining  the  expected 
and/or  high  probability  values  of  U(P±)  for  random  P±. 

It  is  not  difficult  to  show  that  the  expected  value  of  0(P±)  is  Simply  consider 

the  pluses  and  minuses  that  lie  in  the  upper  half  of  the  square.  Since  pluses  in  the  upper 
half  of  the  square  can  only  be  matched  to  minuses  in  the  upper  half  of  the  square,  the 
number  of  upper-half  pluses  less  the  number  of  upper  half  minuses  is  always  a  lower  bound 
on  the  number  of  unmatched  pluses  in  a  maximum  up-right  matching.  The  numbers  of 
upper-half  pluses  and  minuses  are  governed  by  simple  binomial  distributions  with  mean 
N/2  and  variance  N/4.  Hence,  the  number  of  upper-half  pluses  will  exceed  the  number  of 
upper-half  minuses  by  at  least  Cl(y/N)  with  some  constant  positive  probability.  Thus,  the 
expected  value  of  U{P±)  is  at  least  fi(v/N). 

Karp,  Luby  and  Marchetti-Spaccamela  [10]  were  the  first  to  explicitly  define  the  max¬ 
imum  up-right  matching  problem  (i.e.,  the  problem  of  determining  the  expected  value  of 
U(P±)  ).  By  applying  Hall’s  Theorem  [9]  and  analyzing  the  dual  problem,  they  showed 
that  U(P*)  <  0(s/N  log  N)  with  high  probability.  By  applying  an  elegant  argument  of 
Ajtai,  Koml6s  and  Tusnddy  [1],  they  also  showed  that  U(P±)  >  fl(\/jVVlog  N)  with  high 
probability. 

The  resulting  0(>/log  N)  gap  in  the  bounds  remained  until  Shor  [21,  22]  proved  that 
U(P±)  i  fl (y/N  logS/,<  N)  with  very  high  probability.  Not  surprisingly,  this  result  combined 
with  the  0(y/N  log  N)  upper  bound  1>  d  to  the  conjecture  that  U(P±)  =  Q(\/N  logN) 
with  high  probability.  As  it  turns  out,  however,  this  is  not  the  case.  Instead,  it  is  the 
n ( y/N  log3'4  N )  lower  bound  that  is  correct.  In  fact,  we  prove  in  this  paper  that  [/(P*)  = 
0  ( >/iV  log8/4  N)  with  very  high  probability. 

Judging  from  the  similarity  between  the  claimed  bounds  for  the  minimax  grid  matching 
problem  and  the  maximum  up-right  matching  problem,  it  should  no.  be  surprising  that 
the  two  problems  are  closely  related.  In  fact,  any  very  high  probability  upper  bound  on 
L{P)  can  be  transformed  into  a  very  high  probability  upper  bound  for  by  simply 

multiplying  by  0(VN).  To  see  this,  consider  an  up-right  matching  problem  with  N  random 
pluses  P+  and  N  random  minuses  P~ .  (Here  P*  =  P+  U  P~.)  Let  d(N)  be  a  very  high 


probability  upper  bound  on  L(P).  Then,  with  very  high  probability,  L(P+)  <  d(N)  and 
L(P~)  <  d(N).  In  other  words,  the  pluses  and  minuses  can  each  be  matched  to  N  regularly 
spaced  grid  points  so  that  each  grid  point  is  matched  to  a  plus  and  a  minus  that  are  within 
distance  at  most  d(N).  Next  form  a  matching  on  P±  by  matching  the  plus  identified 
with  grid  point  (*,  j)  to  the  minus  identified  with  grid  point  (i  +  2 d(N),j  +  2d(N))  for 
{(*,  j)\l  <  i  <  \/N  —  2 d(N)t  1  <  j  <  y/N  —  2 d(N)}.  For  example,  see  Figure  5. 


Figure  5:  The  matching  of  a  plus  to  a  minus  via  association  with  grid  points.  Since  the 
solid  matching  lines  have  distance  at  most  d(N),  the  minus  is  above  and  to  the  right  of  the 
plus. 

From  inspection  of  Figure  5,  it  is  clear  that  the  procedure  just  described  forms  an  up¬ 
right  matching.  The  only  pluses  not  matched  by  the  procedure  are  those  identified  with 
grid  points  in  the  topmost  2 d(N)  rows  and  the  rightmost  2 d(N)  columns.  As  there  are  less 
than  4d(N)y/N  such  points,  we  can  conclude  that  U(P±)  <  0(y/Nd(N))  with  very  high 
probability.  Hence  the  claimed  very  high  probability  0(y/N  log*/4  N)  bound  for  U{P±) 
immediately  follows  frcn  the  very  high  probability  O(logs^4  N)  bound  for  L(P)  proved  in 
Section  4. 

It  is  conceivable  that  a  symmetric  condition  is  also  true:  namely,  that  a  very  high 
probability  upper  bound  for  U(P *)  can  be  easily  transformed  into  a  very  high  probability 
upper  bound  for  L(P).  We  doubt  that  this  is  the  case,  however,  since  the  minimax  grid 
matching  problem  appears  (at  least  to  us)  to  be  fundamentally  harder  and  more  general 
than  the  maximum  up-right  matching  problem*.  Nevertheless,  the  upper  bound  for 
U(P±)  obtained  from  L(P)  is  tight  up  to  constant  factors. 

2.2.  Two-Dimensional  Discrepancy  Problems 

Discrepancy  problems  have  a  long  and  rich  history  in  mathematics  and  statistics  [6,  19, 
23,  25).  In  fact,  much  of  probability  theory  and  mathematical  statistics  involves  estimating 

*  This  is  true  in  higher  di  .-><*>  ions.  In  d-dimensional  space,  for  d  >  3,  Karp  et  al.  [10]  prove  that  for 
up-right  matching,  the  log  factor  disappears  and  U(P±)  =  8 (N1~xld)\  whereas  there  is  an  easy  lower 
bound  of  nflog^'AO  for  d-dimensional  minimax  grid  matching.  Hence  the  relationship  between  V(Pi) 
and  L(P)  breaks  down  in  higher  dimensions. 


the  likelihood  of  certain  discrepancies  (or  deviations  from  the  norm)  occurring  in  random 
data.  As  a  familiar  example,  consider  an  unbiased  random  walk  of  length  N.  It  is  a  well- 
known  consequence  of  the  reflection  principal  that  the  largest  deviation  from  the  origin 
during  the  walk  is  expected  to  be  Q(\/N).  Hence,  the  largest  “swing”  observed  in  any 
subwalk  of  a  walk  of  length  N  is  also  expected  to  be  0(\/N). 

This  basic  result  has  many  consequences,  and  can  be  stated  more  generally  as  an 
example  of  a  1-dimensional  discrepancy  problem.  For  example,  consider  a  fixed  interval  of 
length  N  that  contains  N  random  points  selected  according  to  independent  and  identical 
uniform  distributions.  Define  the  discrepancy  of  any  subinterval  to  be  the  absolute  value 
of  the  difference  between  the  expected  number  of  points  contained  in  the  interval  (i.e.,  the 
length  of  the  interval)  and  the  actual  number  of  points  contained  in  the  interval.  Using 
the  reflection  principle,  it  is  not  difficult  to  show  that  the  expected  maximum  discrepancy 
over  all  subintervals  is  0(\/N).  Moreover,  with  high  probability  every  subinterval  has 
discrepancy  at  most  0{\[L\J log  N  +  logiV),  where  L  is  the  length  of  the  subinterval. 

In  this  paper,  we  are  interested  in  two-dimensional  generalizations  of  the  preceding 
one-dimensional  results  on  intervals.  Several  different  generalizations  have  been  studied, 
depending  on  how  the  subintervals  are  generalized.  For  example,  consider  a  square  with 
area  N  that  contains  N  random  points  generated  according  to  independent  and  identical 
uniform  distributions.  For  any  region  ft  of  the  square,  define  the  discrepancy  of  the  region 
A  (ft)  to  be  the  absolute  value  of  the  difference  between  the  expected  number  of  points 
contained  in  ft  (i.e.,  the  area  of  ft)  and  the  actual  number  of  points  in  ft.  In  a  classic  paper, 
Kiefer  [11]  showed  that  the  expected  maximum  discrepancy  of  any  oriented  rectangle  is 
0(v^V)  for  fixed  N.  Generalizing  the  law  of  the  iterated  logarithm,  he  also  showed  that 
the  limsup  of  the  maximum  discrepancy  of  any  oriented  rectangle  over  all  N  — ►  oo  is 
0(V^v/log  log  N)  with  probability  1  —  o(i).  These  results  were  later  extended  by  Philipp 
[16]  who  proved  the  same  bounds  for  arbitrary  convex  regions,  and  by  Leighton  [12]  who 
showed  that  with  high  probability,  every  convex  region  has  discrepancy  0{\/A\f\ogN  + 
log  N)  where  A  is  the  area  of  the  region. 

Of  greatest  interest  in  this  paper,  however,  is  the  work  on  maximum  discrepancies  of 
up-right  regions  (also  known  as  lower  layers)  in  the  square.  Up-right  regions  are  defined 
by  monotone  decreasing  curves.  For  example,  Figure  6  illustrates  a  rectilinear  up-right 
region  with  N  =  8,  area  2  and  discrepancy  I. 

By  applying  Hall’s  Matching  Theorem  in  the  usual  way,  it  is  not  difficult  to  see  that  the 
up-right  discrepancy  problem  is  closely  related  to  the  up-right  matching  problem.  More 
precLely,  if  we  consider  a  square  with  area  N  that  contains  N  random  pluses  P+  and  N 
random  minuses  P~ ,  and  we  define  A1  (ft)  to  be  the  number  of  pluses  in  ft  less  the  number 
of  minuses  in  ft,  then  the  maximum  value  of  A*  (ft)  over  all  up-right  regions  is  precisely 
the  number  of  unmatched  pluses  in  a  maximum  up-right  matching  for  P±  —  P+  U  P~ . 
For  random  distributions,  the  maximum  values  of  A* (ft)  and  A(ft)  are  within  constant 
factors  of  each  other  with  very  high  probability,  so  we  can  conclude  from  the  bounds  on 
U(P£)  cited  in  Section  1  that  the  maximum  value  of  A(ft)  over  all  up-right  regions  is 
0(v//Vlog3/4  N)  with  very  high  probability. 

Asymptotically  characterizing  the  expected  maximum  value  of  A  (ft)  over  all  up-right 
regions  has  been  an  open  problem  in  mathematical  statistics  for  some  time.  Perhaps  Blum 
[3]  was  the  first  to  study  the  problem,  proving  a  o(N)  upper  bound  in  1955.  This  was  later 
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Figure  6:  An  up-right  region  with  N  —  8,  area  2  and  discrepancy  1. 

improved  to  0(N3^4)  by  Steele  (23]  in  1977  and  then  to  0(y/N  log  N(\og  log  N)a)  for  some 
constant  a  by  Philipp  [17].  The  first  nontrivial  lower  bound  was  proved  by  Dudley  [7]  in 
1982.  He  showed  that  the  expected  maximum  value  of  A (R)  for  all  up-right  regions  is  at 
least  n (y/N y/\ogN/(\og  log  N)  a +<)  for  any  constant  e  >  0.  Karp,  Luby  and  Marchetti- 
Spaccamela  [10]  were  the  first  to  bring  the  problem  into  the  computer  science  community, 
proving  an  upper  bound  of  0(\Tn  log  N)  and  (using  an  elegant  argument  of  Ajtai,  Komlos 
and  Tusn&dy  [1])  a  lower  bound  of  Cl(\/N y/log  N).  The  0(y/NlogN)  upper  bound  also 
follows  directly  from  prior  work  of  Leighton  and  Leiserson  [13]  and  Ajtai,  Komlds  and 
Tusn&dy  [l].  The  lower  bound  was  later  improved  to  Cl(y/N  \og3^4  N)  by  Shor  (21,  22], 
who  independently  developed  a  method  similar  to  that  followed  by  Dudley.  The  results  in 
this  paper  provide  final  resolution  to  the  problem  by  establishing  a  very  high  probability 
0(\/N  log3'4  N)  upper  bound  on  the  maximum  discrepany  of  any  up-r^ht  region. 

A  priori,  it  is  not  clear  why  we  should  restrict  ourselves  to  either  convex  or  up-right 
regions  of  the  square  when  generalizing  the  1-dimensional  discrepancy  problem.  For  ex¬ 
ample,  why  not  consider  any  simply  connected  region?  (A  region  is  said  to  be  simply 
connected  if  it  is  connected  and  contains  no  holes.)  The  answer,  of  course,  is  that  we  can 
always  construct  a  simply  connected  region  that  contains  all  N  points,  but  that  has  area 
arbitrarily  close  to  zero.  Hence,  the  worst  possible  discrepancy  is  always  achieved. 

Not  all  is  lost,  however,  since  the  discrepancy  of  a  simply  connected  region  can  be 
bounded  in  terms  of  its  perimeter.  In  fact,  we  will  prove  in  Section  4  that  with  very  high 
probability,  the  discrepancy  of  every  simply  connected  region  is  at  most  0(p  log3'4  N  + 
log3/2  N),  where  p  is  the  perimeter  of  the  region.  Note  that  this  result  naturally  generalizes 
the  bound  on  discrepancies  of  up-right  regions  since  p  =  0(y/N)  for  any  up-right  region. 
More  importantly,  however,  the  resuit  will  be  sufficient  to  prove  the  0(log3/4  N)  upper 
bound  for  minimax  grid  matching.  Thu  proof  of  this  fact  is  not  difficult  and  again  uses 
Hall’s  Theorem.  The  details  are  deferred  until  Section  4. 

Although  the  topic  is  not  of  direct  concern  in  this  paper,  two-dimensional  discrepancy 
problems  have  also  been  studied  in  a  worst-case  setting.  The  most  notable  result  in  this 
area  is  due  to  Schmidt  [20]  who  resolved  a  decades  old  open  question  by  proving  that 


no  matter  how  N  points  are  arranged  in  a  square  of  area  N,  there  is  always  an  oriented 
rectangle  in  the  square  which  contains  ft(log.N)  fewer  or  ft  (log  N)  more  points  than  its 
area.  In  other  words,  there  is  always  an  oriented  rectangle  with  discrepancy  ft(log  N). 
Constructions  which  achieve  a  maximum  discrepancy  of  O(logN)  for  this  problem  were 
known  long  ago  (14,  24}  and  have  recently  been  rediscovered  in  the  computer  science 
literature  (8,  15].  Curiously,  it  is  still  not  known  how  to  prove  either  the  “ft(logiV)  fewer” 
or  the  “ft(logAr)  more”  result  individually.  Even  more  startling  is  the  fact  that  there  i 
always  an  oriented  right  triangle  in  the  square  with  discrepancy  ft(iV1/4)  [19]! 

3.  Applications 

In  this  section,  we  briefly  describe  the  applications  of  our  work  to  problems  involving 
bin  packing,  dynamic  allocation,  wafer-scale  integration,  and  testing  pseudorandom  num¬ 
ber  generators.  The  applications  range  widely  in  difficulty,  and  some  are  quite  elegant 
Due  to  space  limitations,  however,  we  will  mostly  just  refer  the  reader  to  relevant  papers 
in  the  literature. 

3.1.  Wafer-Scale  Integration 

The  minimax  grid  matching  problem  was  first  studied  in  the  context  of  wafer-scale 
integration  of  two-dimensional  systolic  arrays  [13],  Only  later  was  its  association  with 
up-right  matching  noticed. 

When  constructing  two-dimensional  arrays  on  a  single  wafer,  the  designer  is  typically 
presented  with  a  regular  y/N  x  y/N  array  of  chip-size  cells,  some  fraction  of  which  are 
functional.  The  task  is  to  connect  the  functional  cells  into  a  smaller  square  array  in  a  way 
that  minimizes  the  length  of  the  longest  wire  needed  to  connect  adjacent  cells  in  the  array. 
For  example,  Figure  7  illustrates  the  connection  of  a  3  x  3  array  of  functional  cells  from  a 
4x4  array  containing  9  functional  cells. 


Defective 

cell 


Figure  7:  Connection  of  a  3x3  array  of  functional  cells. 

Unfortunately,  the  designer’s  task  is  NP-hard  in  general.  However,  the  impact  of  this 
difficulty  is  mitigated  by  two  practical  considerations.  First,  the  defective  cells  tend  to 
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be  located  randomly.  Second,  an  algorithm  need  only  work  well  for  most  (not  all)  wafers. 
Hence,  an  algorithm  which  works  with  high  probability  for  randomly  located  faults  is  very 
acceptable  in  practice. 

The  minimax  grid  matching  problem  suggests  exactly  such  an  algorithm.  First  match 
the  functioning  cells  one-to-one  with  an  imaginary  regularly  spaced  grid  of  the  same  car¬ 
dinality  in  a  way  that  minimizes  the  longest  matching  length.  By  the  work  in  this  paper, 
this  length  will  be  0( log3/4  N)  with  very  high  probability. 

The  one-to-one  association  of  functioning  cells  with  imaginary  grid  points  assigns  a 
label  (*,j)  to  each  functioning  cell.  It  remains  only  to  hook  up  cells  whose  labels  differ 
by  one  in  precisely  one  coordinate.  Since  the  distance  between  any  two  adjacent  cells  (say 
(i,  j)  and  (i,j  +  l),  for  example)  is  at  most  the  distance  between  (*,  j)  and  its  matching  grid 
point  (this  is  0(log3/4  JV)  with  very  high  probability)  plus  the  distance  between  (*, j  +  1) 
and  its  matching  grid  point  (also  0( log8/4  N)  with  very  high  probability)  plus  the  distance 
between  the  two  matching  grid  points  (constant),  every  wire  will  have  length  0(log3/4  N) 
with  very  high  probability. 

Of  course,  this  analysis  assumes  nearly  widthless  wires,  but  for  some  applications  this 
is  reasonable.  Previously,  the  best  upper  bound  proved  for  this  problem  was  O(logJV)  by 
Leighton  and  Leiserson  [13].  The  best  lower  bound  known  remains  the  trivial  bound  of 

n(v/l^lvf). 

3.2.  Two-Dimensional  Bin  Packing 

The  maximum  up-right  matching  problem  was  first  defined  in  the  context  of  two- 
dimensional  bin  packing  [10] .  In  fact,  the  algorithm  proposed  by  Karp,  Luby  and  Marchetti- 
Spaccamela  for  two-dimensional  bin  packing  specifically  uses  the  maximum  up-right  match¬ 
ing  algorithm  as  a  subroutine.  In  what  follows,  we  briefly  describe  the  problem  and  the 
results.  We  refer  the  interested  reader  to  [10]  for  a  description  of  the  algorithm  and  the 
details  of  the  analysis. 

Given  a  collection  of  N  items  from  [0,1]  x  [0,1],  the  two-dimensional  bin  packing 
problem  is  to  pack  the  N  items  into  the  minimum  possible  number  of  unit-square  bins. 
Of  course,  the  two-dimensional  bin  packing  problem  is  AfP-complete.  However,  when 
the  N  items  are  chosen  independently  and  uniformly  from  [0,1]  x  [0,1],  the  Karp-Luby- 
Marchetti-Spaccamela  algorithm  uses  j  -f  0(W)  bins  and  wastes  0(IP)  space  with  very 
high  probability,  where  W  =  Q(y/N  log3/4  N)  is  the  anticipated  number  of  unmatched 
pluses  in  a  random  maximum  up-right  matching  problem.  This  compares  favorably  with 
the  U(y/N)  lower  bound  on  expected  wasted  space  for  any  two-dimensional  bin  packing 
algorithm.  Whether  or  not  there  is  an  algorithm  which  achieves  o(y/N  log3/4  N)  expected 
wasted  space  remains  unknown. 

3.3.  One-Dimensional  Bin  Packing 

Somewhat  surprisingly,  the  apparently  two-dimensional  up-right  matching  problem 
also  arises  in  the  context  of  one-dimensional  bin  packing.  In  particular,  Shor  [21,  22]  showed 
that  the  space  wasted  by  the  Rest  Fit  bin  packing  heuristic  when  packing  N  items  uniformly 
selected  from  [0,l]  is  O(H')  with  high  probability,  where  again  W  =  0(>/Af  log3/4  N)  is 


the  anticipated  number  of  unmatched  pluses  in  a  random  maximum  up-right  matching 
problem.  Moreover,  the  same  upper  bound  holds  when  the  items  are  selected  from  any 
distribution  on  {0,1}  that  is  symmetric  about  1/2. 

The  Best  Fit  algorithm  is  one  of  the  most  common  on-line  algorithms  used  in  practice. 
The  algorithm  packs  each  item  as  it  "arrives”  (never  looking  ahead  at  subsequent  items)  in 
the  fullest  bin  into  which  it  fits.  Although  better  off-line  algorithms  are  known  (First  Fit 
Decreasing  achieves  the  optimal  bound  of  Q(y/N)  expected  wasted  space  [2]),  no  better  on¬ 
line  algorithm  is  known.  In  fact,  Shor  (21,  22]  used  a  related  matching  problem  to  prove 
an  n(v/AT \/log  N)  lower  bound  on  the  expected  wasted  space  of  any  on-line  algorithm, 
assuming  that  the  number  of  items  N  is  not  known  in  advance. 

3.4  Dynamic  Allocation 

The  up-right  matching  problem  also  arises  in  the  context  of  dynamic  algorithms  for 
file  storage.  In  particular,  the  Coffman- Leighton  [4j  Best-Fit-Aligned  (BFA)  algorithm  for 
dynamic  allocation  wastes  0(W)  space  with  very  high  probability  where  W  is  the  expected 
number  of  unmatched  pluses  in  an  N- point  random  up-right  matching  problem,  when  files 
arrive  and  depart  according  to  a  Poisson  process  where  N  is  the  expected  number  of  files 
in  storage  at  any  point  in  time.  Since  the  amount  of  used  (i.e.,  occupied)  space  in  storage 
at  any  time  is  easily  bounded,  this  means  that  the  capacity  of  the  storage  device  can  be 
predetermined  so  as  to  almost  never  waste  more  than  O(v^7log3^4  N)  space  and  so  as  to 
almost  never  need  compaction.  Moreover,  the  results  hold  for  any  distribution  of  file  sizes 
in  (0,1).  These  results  are  much  better  than  those  for  the  worst-case  algorithms  used  in 
practice,  which  commonly  waste  B(N)  space  and  spend  a  constant  fraction  of  their  time 
in  amortized  compaction. 

Although,  we  do  not  have  room  to  describe  the  BFA  algorithm  in  detail  here,  we  can 
mention  a  simpler  balls  and  boxes  problem  with  related  performance.  Consider  an  infinite 
collection  of  boxes  in  a  line  1,2,3,....  Assume  initially  that  precisely  the  first  N  boxes 
are  filled.  At  each  step,  remove  one  of  the  balls  at  random,  and  then  try  to  insert  a  new 
ball  into  one  of  the  first  N  boxes  (chosen  at  random).  Boxes  can  contain  at  most  one  ball 
each,  so  if  we  are  unable  to  insert  the  new  ball  into  its  desired  box,  then  insert  the  ball  into 
the  leftmost  empty  box  to  the  right  of  the  originally  desired  box.  Continue  this  process 
indefinitely. 

It  is  clear  that  the  system  just  described  always  contains  exactly  N  balls,  but  the 
position  of  the  rightmost  ball  will  vary  with  time.  The  question  of  interest  in  dynamic 
allocation  involves  the  expected  location  of  the  rightmost  ball.  This  is,  of  course,  identical 
to  N  plus  the  number  of  empty  boxes  before  the  rightmost  ball  (i.e.,  the  wasted  space). 
By  a  rather  complicated  reduction  to  maximum  up-right  matching,  Coffman  and  Leighton 
|4]  showed  that  the  rightmost  position  at  any  time  T  >  N  log  N  is  N  +  N) 

with  very  high  probability.  Without  up-right  matching  to  rely  upon,  this  result  might  well 
have  been  impossible  to  obtain. 

3.5  Testing  Pseudorandom  Number  Generators 

All  of  the  applications  discussed  thus  far  are  directly  tied  to  one  of  the  matching 
problems.  The  corresponding  discrepancy  problems  also  have  applications,  but  most  are 
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to  problems  in  mathematical  statistics  (e.g.,  see  [23,  26])  and  are  of  limited  interest  to 
computer  scientists.  There  is  one  potential  application  that  could  be  of  interest,  however: 
using  the  minimax  grid  matching  problem  to  test  pseudorandom  number  generators.  As 
an  example,  consider  the  plot  of  “random  numbers”  from  the  IBM  PC  random  number 
generator  shown  in  Figure  8. 


Figure  8:  Plot  of  £000  random  pairs  (Xf,£t+i)>  (x<+4>x>+s)>  •••  generated  from  the  IBM 
PC  Random  Number  Generator  with  IBM  PC  Basie.  (Taken  from  [18].) 

The  points  in  Figure  8  (which  is  copied  from  [18])  were  obtained  by  taking  consecutive 
pairs  of  “random  numbers”  and  using  them  as  z  and  y  coordinates.  In  this  case,  it  is 
readily  apparent  that  the  data  are  not  random.  One  way  of  detecting  this  formally  is 
to  observe  that  there  are  simply  connected  regions  with  area  A  and  perimeter  p  that 
contain  far  more  than  A  +  0(p  log3/4  N)  points.  (Just  draw  the  boundary  around  one  of 
thj  diagonal  clusters.)  Hence,  by  the  arguments  of  Sections  2.2  and  4,  the  minimax  grid 
matching  length  for  this  set  of  points  would  be  far  in  excess  of  0(log3^4  N).  Hence,  we 
would  have  formal  grounds  for  deciding  that  the  IBM  PC  random  number  generator  is 
faulty. 

There  are  several  good  reasons  to  use  the  minimax  grid  matching  length  as  a  barom¬ 
eter  for  randomness.  First,  it  is  a  polynomial  time  test  which  in  some  sense  is  doing  an 
exponential  amount  of  work.  Because  of  Hall’s  Theorem,  we  know  that  if  the  minimax 
matching  length  is  d,  then  every  simply  connected  region  in  the  square  (of  which  there  are 
more  than  an  exponential  number)  is  verified  to  contain  about  the  right  number  of  random 
points  (i.e.,  between  A  —  0(pd)  and  A  +  Q(pd)).  Hence  large  discrepancies  in  odd  shaped 
regions  will  be  readily  detected,  even  if  they  are  not  easily  characterizable  or  readily  appar¬ 
ent  to  the  naked  eye.  Second,  the  test  will  also  identify  generators  that  are  “too  regular.” 
If  the  numbers  are  too  regular,  they  will  form  a  nearly  perfect  grid  resulting  in  a  minimax 
matching  length  that  is  too  small.  Third,  the  distribution  of  the  minimax  matching  length 
is  very  sharply  peaked  at  0(log*^4iV).  The  probability  of  deviating  from  the  peak  by  even 
a  constant  factor  is  vanishingly  (more  than  polynomially)  small.  Hence  small  deviations 
from  the  peak  would  provide  high  confidence  that  the  pseudorandom  number  generator  is 
not  random. 
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The  major  drawback  to  using  the  minimax  matching  length  as  a  statistical  test  is  that 
we  have  not  been  able  to  completely  characterize  its  distribution.  This  would  have  to  be 
done  in  order  for  the  test  to  be  usable  for  fixed  values  of  N.  In  the  mean  time,  the  test 
can  only  be  used  as  a  well-motivated  heuristic.  In  fact,  we  are  currently  doing  just  that. 
Some  students  at  MIT  are  in  the  process  of  coding  several  tests  and  pseudorandom  number 
generators  of  various  types  to  see  how  they  compare. 

The  idea  of  using  a  matching  problem  as  a  test  for  a  pseudorandom  number  generator 
is  not  completely  new.  The  same  idea  (but  for  a  different  problem)  motivated  the  work 
of  Ajtai,  Komlds  and  Tusnddy  (l).  Unfortunately,  they  also  were  unable  to  completely 
characterize  the  distribution  around  the  peak.  More  generally,  related  methods  can  be 
found  in  the  literature  on  spectral  and  Fourier  tests  [5]. 


4.  Proofs 

In  this  section,  we  prove  that  the  minimax  grid  matching  length  L{P)  for  N  random 
points  P  is  O(log3/4  N)  with  very  high  probability.  The  section  is  divided  into  four  subsec¬ 
tions.  In  Section  4.1,  we  formally  convert  the  minimax  grid  matching  problem  into  a  dual 
discrepancy  problem  with  a  straightforward  application  of  Hall’s  Theorem.  We  prove  the 
necessary  bounds  for  the  discrepancy  problem  in  Section  4.3.  Section  4.2  provides  some 
intuition  and  motivation  for  the  rather  complicated  proof  in  Section  4.3,  and  gives  some 
insight  into  why  the  eventual  answer  is  ©(log8/4  JV)  instead  of  ©(^/log  N)  or  G(logW).  In 
Section  4.4,  we  extend  the  result  of  Section  4.3  to  arbitrary  regions  in  the  plane.  The 
result  is  not  needed  for  the  rest  of  the  paper,  but  is  mathematically  more  interesting  and 
natural. 


4.1  Conversion  to  the  Dual  Discrepancy  Problem 

In  the  next  few  subsections,  we  will  restrict  our  attention  to  simply  connected  regions 
with  a  special  rectilinear  form.  In  particular,  we  define  a  partition  T  of  the  square  with 
area  N  into  subsquares,  each  with  side  length  log8/4  N.  In  Section  4.3,  we  will  prove 

that  there  is  a  constant  c  >  0  such  that  with  very  high  probability,  the  discrepancy  of  every 
simply  connected  region  whose  boundary  lies  along  the  edges  of  T  is  at  most  cp  log8/4  N 
where  p  is  the  perimeter  of  the  region.  As  a  consequence,  we  can  conclude  that  with  very 
high  probability,  the  discrepancy  of  every  region  (not  necessarily  connected  or  simple) 
whose  boundary  lies  along  the  edges  of  F  is  at  most  cplog8/4  N  where  p  is  the  perimeter 
of  the  region.  (This  corollary  is  easily  proved  by  decomposing  an  arbitrary  region  into  the 
sum  and  difference  of  simply  connected  regions.)  In  this  subsection,  we  show  how  to  use 
this  result  to  prove  that  the  minimax  grid  matching  length  L(P)  for  N  random  points  P 
is  at  rno3t  d  —  0( log8/4  N)  with  very  high  probability. 

By  Hall’s  Theorem  [9],  it  is  sufficient  to  show  that  for  every  set  A  C  P  with  z  random 
points,  there  are  at  least  z  grid  points  within  distance  d  of  A.  To  do  this,  we  define  a 
coarser  partition  T'  of  the  square  into  ^  subsquares  with  side  length  |  where  (for  now) 
d  —  8c  log3/4  N.  For  any  subset  A  C  P,  let  the  region  R  consist  of  all  the  subsquares  in  V 
containing  a  point  of  A.  Let  R1  be  the  slightly  larger  region  formed  from  R  by  adding  an 
isosceles  right  triangle  with  hypotenuse  ^  to  every  —length  segment  along  the  perimeter  of 
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R.  For  example,  see  Figure  9.  (For  now,  we  overlook  the  special  case  when  the  boundary 
of  R  coincides  with  the  boundary  of  the  overall  square.) 


Added 

triangle 


Figure  9:  The  extension  of  R  to  R1. 

By  construction,  every  grid  point  of  R'  is  within  distance  d  of  some  point  in  A.  Hence 
it  remains  only  to  show  that  the  number  of  grid  points  in  R1  is  at  least  as  large  as  the 
number  of  random  points  in  A.  This  can  be  accomplished  by  observing  that  the  number 
of  grid  points  in  R1  is  at  least  the  area  of  R1  which  is 


Area(fi)  + 


Per(R) 


whereas 


=  Area[R)  +  Per(J2)clogs/4  N 

\A\  <  Area(R)  +  A(P) 

<  Area(R )  +  Per(R)c  log8/4  N  . 


This  completes  the  proof.  (Accounting  for  the  case  when  the  boundary  of  R  coincides 
with  the  boundary  of  the  square  only  multiplies  d  by  a  constant  factor.) 

4.2  Decomposition  Into  Triangles  —  The  Intuition 

In  Section  4.3,  we  will  prove  that  with  very  high  probability,  every  simply  connected 
region  whose  boundary  lies  along  the  edges  of  T  has  discrepancy  at  most  0(p  log5/4  N) 
where  p  is  the  perimeter  of  the  region.  In  what  follows  we  provide  some  intuition  for  the 
result  by  proving  a  much  simpler  result  with  the  same  log3/4  n  bound.  This  simpler  result 
deals  with  the  decomposition  of  a  polygonal  region  into  triangles  and  is  interesting  in  its 
own  right.  The  result  does  not  imply  the  desired  discrepancy  bound,  but  it  does  give  some 
of  the  basic  ideas  behind  it.  In  Section  4.3,  these  ideas  are  obscured  by  technical  details. 
Theorem  1:  Any  polygonal  region  R  with  n  vertices  and  perimeter  p  can  be  decomposed 
into  a  sum  and  difference  of  triangles  T )  such  that  '£.i\jArea{Tj )  =  0(p log®/4  n). 
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Proof:  The  algorithm  we  use  to  decompose  the  region  into  a  sum  and  difference  of  triangles 
is  simple.  At  each  step,  we  reduce  the  number  of  vertices  of  the  polygon  by  half  by  cutting 
off  n/2  consecutive  triangles.  (See  Figure  10.)  In  particular,  we  cut  off  n/2  triangles 
formed  by  pairs  of  adjacent  sides.  We  concentrate  on  what  happens  when  we  remove  the 
jth  triangle,  which  is  formed,  say,  by  vertices  v<  and  v<+t.  Let  the  distance  between 
t>,_i  and  vi+i  be  dy,  and  let  the  two  edges  and  have  lengths  e,  and  e,+t.  Let 

Cj  =  ti  +  Cj+ 1.  In  removing  the  triangle  we  remove  two  edges  totalling  length 

cy,  and  add  an  edge  of  length  dy.  Thus,  we  reduce  the  perimeter  by  6,  =  c,  —  dy. 


Vi 


Figure  10:  Removing  four  triangles. 

We  must  show  that  this  algorithm  produces  a  decomposition  into  triangles  such  that 
£y  yjAreaiTj)  <  0(p  log3/4  n).  Removing  the  jth  triangle  reduces  the  perimeter  by  6 ,  — 
Cj  —  dj.  The  area  of  this  triangle  is  at  most  jdyy^c*  —  dy ,  since  this  area  is  maximized  when 

ei  ~  ei+l- 

We  therefore  get 


s/Arcafr)  <  Id'/'fc}  - 

=  ~  y)v‘ 

<  cV‘«;/4. 

We  repeat  this  step  of  cutting  off  half  the  vertices  until  we  reduce  the  number  of  sides 
to  2,  and  thus  have  no  area  left.  At  each  step  the  total  sum  of  the  cy’s  is  at  most  the 
perimeter  p,  so  since  there  are  logn  steps,  we  have  cy  <  plogn.  Since  6j  was  the 
change  in  the  perimeter  produced  by  removing  the  jth  triangle,  £y=?  <Sy  <  p.  Thus,  we 
have 

££?  -jAreaiX,)  < 

<  2p  log3/4  n. 

The  second  step  is  a  special  case  of  Holder’s  inequality.  I 
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The  expected  discrepancy  of  a  triangle  T  is  &{\J Arta(T)).  If  we  could  assume  that 
the  average  discrepancy  of  a  triangle  in  the  decomposition  of  the  polygonal  region  R  was 
this  expected  discrepancy,  then  by  Theorem  1  the  total  discrepancy  of  the  region  R  would 
be  0(plogs/4  n).  Proving  this  for  a  general  class  of  regions  R  is  the  most  difficult  part  of 
the  paper  and  is  done  in  Section  4.3. 

The  decomposition  used  in  the  proof  of  Theorem  1  is  basically  the  one  that  we  will  use 
for  the  general  result.  We  will  bound,  with  very  high  probability,  the  total  discrepancy  of 
all  the  n/2*  triangles  removed  at  the  Jbth  stage.  We  do  this  by  counting  the  number  of 
possible  sets  of  triangles  that  could  be  removed  at  the  Jbth  stage. 

Although  we  have  not  investigated  the  matter  formally,  the  bound  in  Theorem  1  is 
almost  surely  tight  in  the  worst  case. 


4.3  Formally  Bounding  the  Discrepancy  —  The  Proof 

In  this  section,  we  prove  the  following  central  result. 

Theorem  2:  Consider  a  set  of  N  points  uniformly  and  independently  distributed  in  the 
y/N  x  y/N  square,  and  define  T  to  be  a  grid  of  squares  with  edge  length  log3/4  N.  Then  there 
is  a  constant  c  such  that  with  probability  at  least  1  —  N~io*l/2  N  t  every  simply  connected 
region  R  whose  boundary  follows  T  has  discrepancy  A(J?)  <  cPer(R)  log3/4  N. 

The  proof  is  divided  into  two  sections,  one  deterministic  and  one  probabilistic.  In  the 
deterministic  section,  we  show  that  the  discrepancy  of  any  R  satisfying  the  hypothesis  can 
be  bounded  by  the  sum  of  the  discrepancies  of  2  log  p  disjoint  regions,  most  of  which  are 
drawn  from  "small  enough”  classes  of  "small  enough”  regions.  In  the  probabilistic  part  of 
the  proof,  we  establish  very  high  probability  bounds  on  the  discrepancies  of  the  regions  in 
the  classes,  thus  obtaining  an  upper  bound  on  the  discrepancy  of  any  R. 


4.3.1  The  Deterministic  Part  of  the  Proof 

The  deterministic  part  of  the  proof  of  Theorem  2  consists  of  showing  that  there  are 
collections  of  regions  Di,,<p  of  the  y/N  x  y/N  square  and  a  constant  Co  independent  of  N 
such  that: 

1)  the  area  of  every  region  in  A,fiP  is  at  most  2~'p3y/s, 

2)  the  number  of  regions  in  At,j>  is  at  most  AT*(alog  JV)r+‘,  and 

3)  for  any  p  and  any  simply  connected  region  R  with  boundary  on  Y  and  perimeter  p  in 
the  y/N  x  y/N  square,  there  exist  real  numbers  Sj, . . . ,  sm,  and  regions  Du . . . ,  Dm 
and  D\, . . . ,  D*m  such  that 

3a)  A,  D-  E  for  1  <  *  <  m, 

3b)  A(i2)<E^,(A(A)  +  A(D;)),and 
3c)  «!  +  •••  +  Sm  <  Co. 

The  9it,iP  are  defined  for  values  of  *  that  are  integers  in  the  range  (l,m),  where  m  = 
log p,  values  of  s  that  are  integer  multiples  of  1/  log  N  in  the  range  [1/  log  N,  AT],  and  values 
of  p  that  are  integer  multiples  of  log*/4  N  in  the  range  [log*/4  N,  N).  In  fact,  we  will  assume 
without  loss  of  generality  that  p  is  a  power-of-two  multiple  of  log3/4  N. 


16 


The  regions  in  a  particular  Vi,t,p  need  not  be  connected  or  simple,  but  they  will  consist 
of  the  sum  and  difference  of  easily  described  polygons.  Very  roughly  speaking,  the  regions 
in  P,  #iP  will  correspond  to  the  union  of  triangles  removed  in  the  (logn  -  t)th  halving 
operation  described  in  Section  4.2,  where  s  is  related  to  the  total  area  of  the  removed 
triangles  and  p  is  the  perimeter  of  the  original  region  being  decomposed.  The  precise 
characterization  of  the  is  fairly  difficult,  and  will  be  described  as  we  go  along. 

The  {Di}  and  {£)?}  for  a  region  R  are  derived  from  a  sequence  of  approximations 
Rm+i  to  R  where  R\  =  0  and  Rm+i  =  R.  To  construct  the  sequence,  we  use  two 
different  sequences  of  piecewise  linear  closed  curves.  The  first  sequence  has  vertices  lying 
on  the  boundary  of  R.  We  call  the  *th  curve  of  this  kind  A,.  The  A,  are  then  approximated 
further  by  closed  curves  Bi  with  vertices  on  a  grid  Gy  If  Bi  is  a  simple  closed  curve,  then 
Ri  is  its  interior.  Although  we  begin  with  a  simple  closed  curve,  neither  the  A,  nor  the 
Bi  approximations  are  necessarily  simple  closed  curves,  which  causes  further  problems, 
forcing  un  to  define  a  region  “enclosed”  by  the  curve  By  This  “enclosed”  region  is  R t, 
which  is  not  necessarily  connected,  even  though  R  is  connected. 

To  obtain  the  approximation  Ay  we  mark  2‘  points  at  equal  distances  along  the  bound¬ 
ary  of  R.  We  call  these  points  a,o,  a,i,  . . .,  a,2,_i.  The  starting  point  a,0  =  a0  is  the  same 
for  all  Ai  and  is  chosen  to  be  a  vertex  of  I\  We  then  join  these  points  in  order  by  edges. 
Half  the  vertices  of  Aj+i  are  also  vertices  of  A,-,  specifically,  ay  =  al+22j.  We  let  the  length 
of  the  edge  between  a,-,,-!  and  ay  be  ey. 

The  curve  B,  is  obtained  by  approximating  the  curve  A*  with  points  on  a  grid  Gy  We 
let  the  jth  vertex  6y  of  Bi  be  the  nearest  grid  point  to  the  jth  vertex  ay  of  A,.  If  several 
grid  points  are  equidistant  to  some  vertex,  we  break  the  tie  arbitrarily. 

The  grid  G<  has  points  spaced  evenly  at  distance  =  p/(2'v1ogn).  The  grid  G,-+l  is 
a  refinement  of  Gy  so  a  fourth  of  the  points  of  G,+  j  are  also  points  of  Gy  We  denote  the 
edge  length  of  G,  by  gt  —  <jf1/2,_1. 

We  have  now  produced  the  curve  Bi.  If  it  is  a  simple  closed  curve,  then  the  area  inside 
it  is  our  ith  approximation  72,-  to  the  region  R.  Otherwise,  we  must  do  some  more  work. 
If  an  area  is  enclosed  by  Bt  twice  or  more  (i.e.,  has  winding  number  >  2),  we  still  wish  to 
count  each  point  inside  it  at  most  once  when  calculating  the  discrepancy.  We  do  this  in 
the  following  manner.  If  the  winding  number  of  a  point  is  positive  with  respect  to  By  we 
include  it  in  our  region.  If  the  winding  number  is  zero  or  negative,  we  do  not  include  it. 
This  gives  the  region  72,  which  we  use  as  our  approximation  to  the  region  R. 

By  the  assumption  that  p  be  a  power-of-two  multiple  of  log*^4  Nt  we  know  that  Am+I 
coincides  with  the  boundary  of  R.  Since  the  vertices  of  T  are  a  subset  of  the  grid  points  in 
Gm+ 1,  this  means  that  Bm+i  is  identical  to  the  boundary  of  R  and  thus  72m+J  =  R.  Since 
Bj  contains  only  2  vertices,  Ri  ■-  0.  Hence,  we  have  constructed  the  desired  sequence  of 
approximations  to  R.  The  {D,}  and  {£>*}  are  now  easily  defined  for  1  <  »  <  m  as  follows: 

Di  =  Ri+i/Ri  and 
A  =  RifRi+i, 


where  the  notation  Si /St  is  used  to  denote  the  set  of  points  contained  in  Si  but  not  in  Sa. 
It  is  readily  observed  that 


A{5j)  <  A(S,)  -1-  A [Si/S2)  +  A(S2/Si) 
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for  any  regions  Si  and  S3.  Thus  for  any  t, 


A  (Ri+i)  -  A  (A)  <  A  (A)  +  A(A*)- 

Summing  over  1  <  *  <  m  and  using  Rm+i  =  R  and  Ri  =  0,  we  find  that 

A(fi)  <  E|A(D.)  +  A(Df)|, 

i=l 

as  desired. 

Each  A  and  D*  is  contained  in  the  union  of  “triangles”  removed  from  B,+i  to  produce 
Bi.  That  is  because  any  point  which  has  a  positive  winding  number  with  respect  to  B,+1  but 
not  with  respect  to  Bi  must  be  contained  in  one  of  the  triangie-like  objects  removed  from 
Bi+i  to  form  Bi.  (The  objects  removed  from  to  form  At  are  triangles,  albeit  possibly 
overlapping,  but  slight  deformations  of  this  basic  structure  are  possible  when  passing  to 
£»+ 1  and  Bi  because  of  the  coarseness  of  the  approximation.)  We  will  eventually  use  this 
fact  to  show  that  the  areas  of  the  A’8  and  D-'s  are  small.  Before  doing  this,  however,  we 
must  classify  A  and  A*  a8  elements  of  Di>t,p  for  some  s  and  p.  The  choice  of  p  is  easy;  it  is 
simply  the  perimeter  of  the  original  region  R.  Since  the  boundary  of  R  follows  F,  p  must 
be  an  integer  multiple  of  log*/4  N  in  the  range  [log*/4  N,  AT),  as  needed  for  the  definition  of 
j  In  fact,  we  can  increase  p  to  a  power-of-two  multiple  of  log*/4  N  by  appending  to  R 
a  zero-area  region  with  sufficiently  large  perimeter. 

The  selection  of  s  =  s,  is  much  more  complicated  and  depends  on  the  triangle  decom¬ 
position  of  R.  In  particular,  we  will  need  to  examine  the  triangles  Tj,  that  are  removed 
from  .Aj+i  to  form  Ai.  Specifically,  we  define  T|;  to  be  the  triangle  with  vertices  a,  ;_i,  atJ 
and  al+i  jy_i.  Recall  that  e,-y  is  defined  to  be  the  distance  between  and  a,  y  and  define 


2‘  * 


j 

•y 


The  value  of  is  a  normalized,  dimensionless  parameter  that  is  closely  related  to  the 
perimeter  of  Ai  as  well  as  the  area  of  i2,+i/A-  It  is  hard  to  provide  better  motivation  for 
this  definition  except  to  say  that  it  has  a  sufficiently  strong  form  so  that  the  inequalities 
will  all  work  out  right  later.  Unfortunately,  the  same  cannot  be  said  of  more  natural 
definitions.  The  following  fact,  in  particular,  will  prove  useful  later. 

Fact  1:  qi  <  qi+i  <  1  for  0  <  *  <  m  —  1. 

Proof:  To  show  that  qm  <  1,  we  observe  that  the  2m  edges  of  Am  each  have  length  at 
most  p/2m.  Thus 

9m  <  — r2m(p/2m)*  =  1. 

P 

To  show  that  <7,  <  g,-+1,  we  apply  the  triangle  inequality  to  TtJ  to  show  that  e*;  < 
2(c?+i,2;-i  +  ei+i,ij)-  Summing  over  1  <  j  <  2'  yields  that 


and  thus  that 


a*  9.+1  a,+l 

V£«‘+I, 

p  i=l  p  i=l 


which  means  that  g,  <  I 

It  is  also  worthwhile  defining 


~  2(ei+l,2j-l  +  Ci+l,lj)  —  e<i 

which  is  proportional  to  the  contribution  of  triangle  Ttj  to  qit  and 


—  9«+l  9i 


which  is  the  difference  in  q  for  A{+\  and  A,-.  By  Fact  1,  r,-  >  0  for  0  <  t  <  m  —  1  and 

m- 1 

Y  fi  ~  9m  <  1. 

i=l 


Also  note  that 


Y  k<i  =  P*r«  /2<- 

i=i 


Lastly,  we  define 


r|  =  r,  + 


ri 

2»‘l 


and  set  s4  to  be  the  smallest  integer  multiple  of  1/  log  N  that  is  greater  than  ar'  +  fi/  log  N 
where  a  and  fi  are  sufficiently  large  constants  to  be  determined  later.  (It  will  turn  out 
that  a  —  9ex  and  fi  =  9e(l  +  16tt)  are  sufficient  where  e  =  2.718....)  The  class  is  then 
defined  to  be  the  set  of  all  Di  and  Df  resulting  from  a  region  R  with  perimeter  p  for  which 
s  =  Si.  It  now  remains  to  prove  the  hypotheses  stated  at  the  beginning  of  this  section.  In 
particular,  we  still  must  show  that: 


1)  Area(Di)  and  Area(D^)  are  at  most  2~,p,v/a4, 

2)  the  number  of  Di  or  D}  for  which  for  some  integer  t  is  at  most 

N2tr+l ,  and 

3)  Si  + - h  sm  <  Co  where  Co  is  a  constant  independent  of  N,  p  and  R. 

We  start  by  showing  that  the  areas  of  Di  and  D,*  are  at  most  2~'ply/si.  From  before, 
we  know  that  Di  and  £>*  are  contained  in  the  union  of  the  triangle-like  objects  removed 
from  Bi+i  to  form  B,-.  In  particular,  define  E{  =  U*Li  Ty  and  define  F<  to  be  the  set  of 
all  points  within  distance  of  a  point  in  B,.  Then  D,-,  Df  C  Fi  and  it  is  sufficient 

to  show  Area(Fi)  <  2"'pJv/s4.  Since  £?,•  can  be  bounded  with  two  closed  curves  each  with 
length  at  most  p,  it  is  easily  seen  that 


Area(Fi)  <  Area{Ei)  +  2\/2flf4+1p  +  2xg]+v. 

The  area  of  Ei  is  bounded  by  the  sum  of  the  areas  of  the  triangles  T4j  it  contains. 
The  area  of  Tijt  in  turn,  is  at  most  56,,  times  the  altitude  of  the  triangle.  A  simple 
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geometric  argument  shows  that  the  altitude  is  at  most  \\fk7j  where,  as  before,  kl}  = 

2(e-+1Jj_1  +  ef+1Jy)  -  efy.  Since  ey  <  p/2\  this  means  that  v4rea(7y)  <  psJk~if2i+t  and 
thus  that 

Area(Ei)  <  - £7  £  ^ 

<  y^u  %/2* 

p\/5 

2<+J  • 

We  can  now  conclude  that 

/trea(Fj)  <  +  2\Z2gi+1  +  2*rft?+1 

<  2-yf—  +  -^-  + - - - ^ 

F  V  4  ploglV  2‘+Mog JV; 

<  2"'p1Y/ii 

provided  only  that  a  >  and  0  >  f .  (Recall  that  ft  =  p/2'v/Iog  IV  and  p  >  log3/4  N  since 
the  boundary  of  R  coincides  with  I\) 

We  next  show  that  the  number  of  D{  and  £><  for  which  <  s<  <  f°r  some  integer 

t  is  at  most  tr+t.  This  will  be  reasoned  with  a  Kolmogorov  complexity  type  of  argument. 
In  particular,  we  will  associate  every  set  of  sequential  approximations  {B,, . . . ,  Bi+1}  for 
which  Sj  <  y  with  a  unique  triple  of  lists.  The  first  list  will  consist  of  two  integers 
from  [1,1V].  The  second  list  will  consist  of  2*+l  integers  drawn  from  [1,9].  The  third  list 
will  consist  of  2<+1  integers  that  sum  to  at  most  ;2,+l(t  —  0  +  16a),  where  a  and  0  are 
the  constants  in  the  definition  of  s<  —  ar{  -f  0/  log  N.  This  will  be  sufficient  to  prove  the 
desired  result  since  the  number  of  such  triples  of  lists  is  at  most 


NW,.^(t-^e«)  +  r" 

=  "£  +  160  +  ■>] 

=  W1[8I!(^-^+16a+;)]1'+, 

<  N't2*' 


provided  that  o  >  9e?r  and  0  >  16a  +  -j-  >  9e(l  +  16x).  Notice  that  we  used  the  well-known 
inequality  <  (7)*  in  the  third  step  of  the  calculation. 
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We  start  describing  the  association  of  {Bi, . . . ,  B,+i}  with  lists  by  considering  the  ways 
that  Bi  can  be  extended  to  form  B<+ 1.  First  we  note  that  for  every  vertex  61}  of  B,,  the 
corresponding  vertex  6,+i^y  of  B*+ i  is  either  the  same  point  or  an  adjacent  point  of  the 
grid  Hence  there  are  9  possible  ways  to  refine  the  approximation  of  each  vertex  of 

B{  in  B,+ 1- 


• 

21  10  14 

20  9  2  6  15 

•••  13  5  1  3  11  ••• 

19  8  4  7  16 

18  12  17 


Figure  11:  Numbering  the  grid  points . 

The  remaining  vertices  of  Bi+l  fall  “between”  consecutive  vertices  of  B<.  To  specify  the 
location  of  vertex  6t+ijy_|,  we  label  ail  the  grid  points  of  Gi+i,  starting  with  the  midpoint 
of  the  edge  6j,y-i6iy  of  B,.  The  labelling  is  done  in  increasing  order  of  distance  from  the 
midpoint.  For  example,  see  Figure  11.  Since  the  edges  of  Ai+l  have  length  at  most  p/2’+1, 
the  next  vertex  of  B,+1  conceivably  could  be  as  far  as  p/2*+1  away  from  this  midpoint  (up  to 
the  0(<7i+1)  error  induced  by  the  approximation).  This  is  essentially  \/log  N  grid  points  of 
Gi+i  away.  However,  most  points  will  be  considerably  closer,  resulting  in  a  much  narrower 
range  of  choices  for  most  6j+i,ay-i.  We  will  prove  this  using  the  following  elementary  fact 
from  geometry. 

Fact  2:  Let  UV  denote  a  line  segment  and  let  \UV\  denote  its  length.  The  locus  of  points 
vV  determined  by  2(\UW\2  +  |VW|J)  —  \UV |2  =  k  is  a  circle  of  radius  |v/fc  about  the 
midpoint  of  UV . 

Idy  Fact  2  the  point  6<+1i2j_i  jf  Bl+j  will  be  in  a  circle  of  radius  p  =  \{\fkij  +  3y/2gi+i) 
centered  on  the  midpoint  of  the  edge  (Recall  that  fcj,  =  2(e?fl  2j_,  +  e*+l  2;)  - 

and  that  B,+j  is  an  approximation  to  Since  we  labelled  the  points  in  order  of 

increasing  distance  from  the  midpoint,  the  label  of  the  point  will  be  less  than 

the  number  of  grid  points  in  that  circle.  The  number  of  grid  points  of  G,+1  in  a  circle  of 
radius  p  is  at  most  i r(p  +  gi+i/\/2)2/g2+l.  It  follows  that  the  label  /,/  of  the  point  6,+i,2>-i 
satisfies  the  equation 


Summing  this  inequality  for  ftJ-  above,  we  get  that 


>=»  i=» 


"  V  2'  I 


+  16jt2' 


=  jr2*+lri  log  N  +  16x2*. 

The  preceding  analysis  indicates  that  the  number  of  options  for  Bi+i  given  B,  is  limited 
in  a  fundamental  way  by  r,-.  In  particular  Bl+ j  can  be  completely  specified  by  B,,  r,, 
a  list  of  2‘  numbers  from  [1,9],  and  a  list  of  2‘  positive  integers  that  sum  to  at  most 
?r2,+1r,  log  JV  -1-  16x2*.  Applying  this  process  recursively,  we  find  that  {B2>. . . ,  Bt+l}  can 


be  completely  specified  by  Bj,  rj,  a  list  of  2*  +  2*-1  -f  . 
and  a  list  of  2*+1  positive  integers  that  sum  to  at  most 


+  2  <  2*+1  numbers  from  [1,9], 


£  (x2*,+1r,«  log  N  +  163t2*#)  <  *r[2,+I  log  N  +  16tt2,+1  . 

•*=i 

If  we  restrict  our  attention  to  those  {B\t . . .  ,  B,-+i}  for  which  a,  <  t/  log  TV,  then  the 
preceding  sum  is  at  most 

x2i+l  ^*log^~  ^  +  16*2<+1  <  ^(t-/?  +  167r). 

Since  Bx  contains  only  two  points  from  Grj,  it  can  be  completely  specified  by  naming 
the  points.  Since  Gx  has  at  most  N  points,  this  is  accomplished  with  a  list  of  two  numbers 
from  [1,2V].  Hence,  every  {Bu . . .  ,B,+1}  for  which  s<  <  can  be  uniquely  specified  by 
a  triple  of  lists:  one  with  two  numbers  from  [1,2V],  one  with  2,+1  numbers  from  [1,9],  and 
one  with  2<+1  integers  that  sura  to  at  most  |2,+l(t  -  0  +  16x).  As  argued  previously,  this 
means  that  the  number  of  B,  and  D-  for  which  a,  <  is  at  most  N2t2'+l . 

The  last  step  in  the  deterministic  part  of  the  proof  is  to  show  that  the  maximum 
of  si  +  •  •  •  +  sm  over  all  regions  R  is  bounded  by  a  constant  independent  of  N  and  p. 
Fortunately,  this  is  easy  to  do  since 
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log  AT 


t 

•=1  \ 


logiV, 


_  (/?  +  !)m  .  y*  y'  ri 

log  n  +  hh 2,w 

=  (/?-H)l°gp  A  A  r, 

log  N  flti  2’W 

m 

<  (3  +  1  +  a  £  2r, 

;= i 

^  3  -f- 1  -J-  2a 


which  is  a  constant  independent  of  N  and  p. 


4.3.2  The  Probabilistic  Part  of  the  Proof 

The  probabilistic  part  of  the  proof  of  Theorem  2  consists  of  bounding  the  discrepancies 
of  the  regions  in  the  P,-#p.  To  do  this,  we  will  make  use  of  the  following  simple  fact,  which 
can  be  proved  by  an  elementary  counting  argument. 

Fact  3:  Given  a  region  D  with  area  A  in  the  y/N  x  y/N  square,  the  probability  that 
A(D)  >  q  is  at  most 

O  for  q  <  A  and 

O  ^  for  q  >  A, 

where  cx  is  a  constant  independent  of  D  and  N. 

In  order  to  bound  the  maximum  discrepancy  of  a  region  in  P,,p  with  very  high  prob¬ 
ability,  it  is  sufficient  to  find  a  q  such  that 

^a^+Mogf.iog  N)e-cxi*/A  <  o(  2-'°«*/aw)  and 

^a22’+ilo*f*lo*;v)e"ei9  <  O  ^2~lot*/awj 

where  A  <  2~‘p2v/s.  This  is  because  the  probability  that  a  fixed  region  in  P,iJp  has 
discrepancy  q  or  greater  is  0(e_e,,3/>4)  or  0(e_•,,)  depending  on  whether  or  not  q  <  A, 
and  because  Pitp  only  contains  at  most  jV222‘t,lo*(*,°*/v)  regions. 

Solving  for  q  =  qx  in  the  first  inequality  we  find  that  it  is  sufficient  for 

qx  =  O  (yjA\\og*l7  N  +  2*  log(slogiV))^ 

<  O  ^2 ~'/Jps!/4  log3/4  N  +  ps1t*yf\og(s  log  iV)^ 

<  O  (2 ~'/2p  log3/4  N  +  ps1/4  log(s  log  N))  . 

Solving  for  q  —  q2  in  the  second  inequality,  we  find  that  it  is  sufficient  for 

q2  =  O(iog‘,‘‘  N  +  2‘  log(slog  N)) 

—  O(!og3/2  N  +  2'  log  log  N). 

Home  with  very  high  probability,  the  maximum  discrepancy  of  a  region  in  Pq,iP  is  at  most 
q  =  9i  +  92 

=  0(2“,/Jp  logs/4  N  +  ps1/4  log(s  log  N)  +  log3/2  N  +  2*'  log  log  N). 

Since  there  are  only  a  polynomial  number  of  combinations  of  i,  s  and  p,  we  can  conclude 
that  the  preceding  bound  holds  with  very  high  probability  for  all  P,,»,p.  From  Section  4.3.1, 
we  know  that  for  every  region  R, 

A(ft)  <  f;(A(D,)+  AfD-)) 

i  = 1 


23 


where  D]  e  for  1  <  i  <  m  and  3i  +  •  •  •  +  sm  <  co  for  some  constant  c0.  Thus, 
with  very  high  probability,  the  discrepancy  of  any  region  R  whose  boundary  follows  T  is 
at  most 


(logp 

53  (2 ~*!%p  log8/4  TV  +  paV*  log(s<  log  TV)  +  log3/1  TV  +  2'  log  log  TV) 
where  si  <  Cq. 

We  will  consider  the  impact  of  each  of  the  four  terms  in  the  sum  individually.  The  first 
term  in  the  sum  forms  a  geometric  series  that  sums  to  0(p log8/4  TV).  The  second  term  in 
the  sum  is  the  most  important  and  lies  at  the  heart  of  the  proof.  Because  x1/4  log(x  log  TV) 
is  convex  in  x,  we  know  by  Jensen’s  inequality  that  the  value  of  pa)**  log (5,  log  TV)  is 
maximized  when  the  s,  are  all  equal.  Hence  these  terms  sum  to 

0(log  TVp(c0/  log  TV)1/4  log  c0)  =  0(p  log8/4  TV) . 

The  third  term  in  the  sum  trivially  adds  up  to  O (log p  log8/1  TV)  which  is  0(p  log8/4  TV) 
provided  that  p  >  log8/4  TV  log  log  TV.  The  fourth  term  in  the  sum  forms  another  geometric 
series  that  converges  to  0(p  log  log  TV). 

The  preceding  argument  concludes  the  proof  of  Theorem  2  except  in  the  special  case 
when  p  <  log8/4  TV  log  log  TV.  For  this  case,  we  can  forget  Section  4.3.1  and  apply  Fact  3 
directly. 

Since  such  regions  can  contain  at  most  (log  log  TV)1  interlinked  Bquare  cells  of  T,  there 
are  at  most  0(TVJ)  such  regions.  Hence,  with  very  high  probability,  the  maximum  discrep¬ 
ancy  of  such  a  region  is  at  most  q  =  q%  +  q2  where 

<  o( 2-|0«*/,Af)  and 

N*e-'in  <  o  (2'l0,*/a*)  . 

Solving  the  equations  yields 

qx  =  0(vC41og8/4  TV) 

<  0(p  log8/4  TV) 

and 

q2  =  O  (log8/1  TV) 

<  0(p  log*/4  TV) 

since  p  >  log3/4  TV  and  A  <  0(p3).  Thus,  with  very  high  probability,  the  discrepancy  of 
any  region  with  p  <  log8/4  log  log  TV  is  0(p  log8/4  TV)  and  we  are  finally  done  with  the  proof 
of  Theorem  2! 

It  is  worth  remarking  that  the  probability  bound  in  Theorem  2  is  essentially  tight. 
For  example,  the  discrepancy  of  the  upper  half  of  the  y/N  x  \/N  square  will  exceed 


ci\/N  log' /4  N  with  probability  2~e'1  logS/3  N  for  some  constant  ci  and  c2.  In  fact,  this  ex¬ 
ample  iihstrates  the  relationship  between  the  constants  in  the  ©(plog3/4  N)  bound  and 
the  constant  in  the  exponent  of  the  very  high  probability  bound.  With  further  work,  it  is 
conceivable  that  the  exact  relationship  could  be  specified. 

4.4  Extension  to  Arbitrary  Regions 

In  what  follows,  we  show  how  to  extend  the  discrepancy  bound  in  Theorem  2  to 
arbitrary  simply  connected  regions.  The  result  is  not  essential  to  the  rest  of  the  paper 
(Theorem  2  is  sufficient),  but  it  is  more  natural  than  Theorem  2  and  may  eventually  prove 
to  be  equally  important. 

Theorem  3:  Consider  a  set  of  N  points  uniformly  and  independently  distributed  in  the 
y7v  x  \fN  square.  There  is  a  constant  c  such  that  with  probability  1  —  N  lo«,/3'v)  every 
simply  connected  region  R  has  discrepancy  A (R)  <  cPer(R)  log3/4  N  +  clog3/2  IV. 

Proof:  The  proof  is  not  difficult  given  Theorem  2  so  we  will  only  sketch  the  argument 
here.  Divide  R  into  two  regions  R'  and  R"  where  R'  consists  of  all  log3/4  N  x  log3/4  N 
square  subregions  of  T  that  are  entirely  contained  in  R  and  where  R"  =  R/R'.  Although 
R’  need  not  be  connected,  we  can  conclude  immediately  from  Theorem  2  that  A(R')  < 
0(p'  log8/1  N)  with  very  high  probability  where  p'  is  the  perimeter  of  R*.  By  elementary 
geometry  p'  <  6p  where  p  is  the  perimeter  of  R  and  hence  A(J2')  <  0(p  log3/4  N)  with 
very  high  probability.  Since  A (R)  <  A(R')  +  A(R"),  it  remains  only  to  show  that  A(R")  < 
0(p  log3/  N  T  log3/2  N)  with  very  high  probability. 

With  very  high  probability,  every  square  subregion  of  T  contains  at  most  0(log3/2  N) 
random  points.  Since  any  square  containing  part  of  R"  must  contain  part  of  the  simple 
closed  cu  :ve  boundary  R,  at  most  0(1  +  p/  log*/4  N)  squares  can  contain  part  of  R".  Any 
such  square  can  contribute  at  most  O(log3/2  N)  to  the  discrepancy  of  R"  since  all  squares 
have  0(log3/2  N)  area  and  0(log3/3iV)  random  points  with  very  high  probability.  Thus 
A  (ft")  <  0(  log3/2  N  +  p  log3/4  IV)  with  very  high  probability,  as  claimed.  I 

It  is  vorth  remarking  that  a  slightly  stronger  result  can  be  proved  if  the  very  high 
probability  assumption  is  relaxed.  In  particular,  with  high  probability,  every  region  with 
pcrimotc-  p  has  discrepancy  at  most  0(log  N  +  plog3/4  N).  (This  can  be  proved  by  noting 
that  Theorem  3  handles  the  case  when  p  >  log3/4  IV,  and  then  observing  that  there  are 
at  most  .i  polynomial  number  of  regions  with  p  <  log3/4  IV.)  The  log3/4  N  factor  cannot 
be  improved  since  Shor  [22]  has  demonstrated  the  existence  of  regions  having  discrepancy 
M(p  log3/'  N)  for  p  =  Q(\fN)  with  very  high  probability. 

5.  Bom. irks 

Sjpoi;  reflection,  it  really  should  not  be  surprising  that  the  minimax  grid  matching 
problem  arises  in  so  many  diverse  and  useful  applications.  After  all,  its  dual  discrepancy 
problem  captures  a  very  important  measure  of  expected  discrepancies  in  random  data.  In 
fart ,  it  i«  precisely  the  measure  that  is  important  to  the  analysis  of  many  algorithms. 

it  is  likely  that  other  applications  of  this  work  will  be  discovered,  and  even  more 
likely  sti'l  that  other  matching  and/or  discrepancy  problems  will  arise  in  the  analysis 
of  algori  hms.  One  such  problem  that  still  remains  unsolved  is  the  rightward  matching 
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problem.  In  rightward  matching,  we  have  N  random  pluses  and  N  random  minuses  P* 
in  a  unit  square,  and  are  asked  to  match  pluses  rightward  to  minuses  in  a  way  which 
minimizes  the  average  vertical  distance  V(P±)  of  the  matching  edges.  (Unmatched  pluses 
are  considered  to  match  to  the  top  or  bottom  of  the  square.)  As  a  consequence  of  the 
minimax  grid  matching  result  in  this  paper,  we  can  conclude  that  V(P±)  <  0(log3^47V) 
with  very  high  probability.  Improvement  of  the  lower  bound  of  Q(yflogN)  would  directly 
lead  to  improved  lower  bounds  on  the  expected  wasted  space  of  any  on-line  bin  packing 
algorithm  [21,  22]. 

Discrepancy  problems  are  still  of  interest  in  mathematical  statistics.  For  example, 
Dudley  [7]  leaves  open  the  question  concerning  the  maximum  expected  discrepancy  of  any 
convex  3-dimensional  region.  Curiously,  this  appears  to  be  analogous  to  the  2-dimensional 
up-right  discrepancy  problem  analyzed  in  this  paper.  We  suspect  that  the  bound  is  again 
Q(y/N  log3/*  N)  but  have  not  proved  it. 

Multidimensional  matching  and  discrepancy  results  will  also  be  likely  to  have  conse¬ 
quences  for  multidimensional  packing  and  allocation  problems.  For  example,  see  [10]. 
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